Palindromic target site identification in SARS-CoV-2, MERS-CoV and SARS-CoV-1 by adopting CRISPR-Cas technique.

Clustered Regularly Interspaced Short Palindromic Repeats (CRISPR) associated Cas protein (CRISPR-Cas) has turned out to be a very important tool for the rapid detection of viruses. This can be used for the identification of the target site in a virus by identifying a 3-6 nt length Protospacer Adjacent Motif (PAM) adjacent to the potential target site, thus motivating us to adopt CRISPR-Cas technique to identify SARS-CoV-2 as well as other members of Coronaviridae family. In this regard, we have developed a fast and effective method using k-mer technique in order to identify the PAM by scanning the whole genome of the respective virus. Subsequently, palindromic sequences adjacent to the PAM locations are identified as the potential target sites. Palindromes are considered in this work as they are known to identify viruses. Once all the palindrome-PAM combinations are identified, PAMs specific for the RNA-guided DNA Cas9/Cas12 endonuclease are identified to bind and cut the target sites. In this regard, PAMs such as 5'-TGG-3' and 5'-TTTA-3' in NSP3 and Exon for SARS-CoV-2, 5'-GGG-3' and 5'-TGG-3' in Exon and NSP2 for MERS-CoV and 5'-AGG-3' and 5'-TTTG-3' in Helicase and NSP3 respectively for SARS-CoV-1 are identified corresponding to SpCas9 and FnCas12a endonucleases. Finally, to recognise the target sites of Coronaviridae family as cleaved by SpCas9 and FnCas12a, complements of the palindromic target regions are designed as primers or guide RNA (gRNA). Therefore, such complementary gRNAs along with respective Cas proteins can be considered in assays for the identification of SARS-CoV-2, MERS-CoV and SARS-CoV-1.

Correction of a Factor VIII genomic inversion with designer-recombinases.

Despite advances in nuclease-based genome editing technologies, correcting human disease-causing genomic inversions remains a challenge. Here, we describe the potential use of a recombinase-based system to correct the 140 kb inversion of the F8 gene frequently found in patients diagnosed with severe Hemophilia A. Employing substrate-linked directed molecular evolution, we develop a coupled heterodimeric recombinase system (RecF8) achieving 30% inversion of the target sequence in human tissue culture cells. Transient RecF8 treatment of endothelial cells, differentiated from patient-derived induced pluripotent stem cells (iPSCs) of a hemophilic donor, results in 12% correction of the inversion and restores Factor VIII mRNA expression. In this work, we present designer-recombinases as an efficient and specific means towards treatment of monogenic diseases caused by large gene inversions.

LIRBase: a comprehensive database of long inverted repeats in eukaryotic genomes.

Small RNAs (sRNAs) constitute a large portion of functional elements in eukaryotic genomes. Long inverted repeats (LIRs) can be transcribed into long hairpin RNAs (hpRNAs), which can further be processed into small interfering RNAs (siRNAs) with vital biological roles. In this study, we systematically identified a total of 6 619 473 LIRs in 424 eukaryotic genomes and developed LIRBase (https://venyao.xyz/lirbase/), a specialized database of LIRs across different eukaryotic genomes aiming to facilitate the annotation and identification of LIRs encoding long hpRNAs and siRNAs. LIRBase houses a comprehensive collection of LIRs identified in a wide range of eukaryotic genomes. In addition, LIRBase not only allows users to browse and search the identified LIRs in any eukaryotic genome(s) of interest available in GenBank, but also provides friendly web functionalities to facilitate users to identify LIRs in user-uploaded sequences, align sRNA sequencing data to LIRs, perform differential expression analysis of LIRs, predict mRNA targets for LIR-derived siRNAs, and visualize the secondary structure of candidate long hpRNAs encoded by LIRs. As demonstrated by two case studies, collectively, LIRBase bears the great utility for systematic investigation and characterization of LIRs and functional exploration of potential roles of LIRs and their derived siRNAs in diverse species.

A Conserved Histophilus somni 23S Intervening Sequence Yields Functional, Fragmented 23S rRNA.

Histophilus somni is a Gram-negative bacterial organism that acts as an opportunistic pathogen and is a fastidious member of the Pasteurellaceae family associated with diseases of respiratory, reproductive, cardiac, and other tissues of ruminants. We identified an intervening sequence (IVS) embedded in all five copies of the 23S rRNA gene in the closed genome sequence of the H. somni isolate USDA-ARS-USMARC-63250 that may play an important role in affecting the biology of the organism. Sequencing the RNA from this isolate shows that most of the IVS is cleaved from the transcript, resulting in independent fragments of this structural rRNA that remain functional within the bacterial ribosome. The IVS lies between positions 1170 and 1278 bp of the 3,017-bp gene and exhibits self-complementarity between its 5' and 3' ends that predicts a stem-loop structure interrupting helix-45 in the transcribed 23S rRNA. Excision removes a 94-nucleotide (nt) stem-loop structure that displays an unusual 1-nt 3' end overhang instead of the more typical 2-nt overhang commonly observed at the ends of other excised IVS stem-loops. A comparison with genomes of other H. somni isolates indicates that this IVS is highly conserved, with 31 of 32 complete genomes having similar interruptions of canonical 23S rRNA genes. The potential biological effects of either the released IVS or the fragmentation of the functional 23S rRNA are unknown, but fragmentation may enhance rRNA degradation in ways that contribute to the regulation of gene expression. IMPORTANCE The genome biology underlying H. somni virulence, pathogenicity, environmental adaptability, and broad tissue tropism is understood poorly. We identified a novel H. somni 109-nt IVS stem-loop structure, of which the central portion is excised from the 23S rRNA transcript, resulting in the fragmentation of this rRNA in the H. somni isolate USDA-ARS-USMARC-63250 and the release of a 94-nt structured RNA of unknown function. We determined that this peculiar rRNA biology is widespread among sequenced H. somni isolates, suggesting it has importance to organism biology. The fragmented 23S rRNA molecules remain functional in the ribosome, given that the isolate grows in culture. The structured excised portion of the IVS, presumably due to the action of the endoribonuclease III, has an unusual 1-nt 3' end overhang. This newly discovered H. somni 23S rRNA fragmentation may enhance rRNA degradation providing a previously unrecognized avenue for regulating H. somni biological processes.

Ligands as Stabilizers of G-Quadruplexes in Non-Coding RNAs.

The non-coding RNAs (ncRNA) are RNA transcripts with different sizes, structures and biological functions that do not encode functional proteins. RNA G-quadruplexes (rG4s) have been found in small and long ncRNAs. The existence of an equilibrium between rG4 and stem-loop structures in ncRNAs and its effect on biological processes remains unexplored. For example, deviation from the stem-loop leads to deregulated mature miRNA levels, demonstrating that miRNA biogenesis can be modulated by ions or small molecules. In light of this, we report several examples of rG4s in certain types of ncRNAs, and the implications of G4 stabilization using small molecules, also known as G4 ligands, in the regulation of gene expression, miRNA biogenesis, and miRNA-mRNA interactions. Until now, different G4 ligands scaffolds were synthesized for these targets. The regulatory role of the above-mentioned rG4s in ncRNAs can be used as novel therapeutic approaches for adjusting miRNA levels.

YB-1 unwinds mRNA secondary structures in vitro and negatively regulates stress granule assembly in HeLa cells.

In the absence of the scanning ribosomes that unwind mRNA coding sequences and 5'UTRs, mRNAs are likely to form secondary structures and intermolecular bridges. Intermolecular base pairing of non polysomal mRNAs is involved in stress granule (SG) assembly when the pool of mRNAs freed from ribosomes increases during cellular stress. Here, we unravel the structural mechanisms by which a major partner of dormant mRNAs, YB-1 (YBX1), unwinds mRNA secondary structures without ATP consumption by using its conserved cold-shock domain to destabilize RNA stem/loops and its unstructured C-terminal domain to secure RNA unwinding. At endogenous levels, YB-1 facilitates SG disassembly during arsenite stress recovery. In addition, overexpression of wild-type YB-1 and to a lesser extent unwinding-defective mutants inhibit SG assembly in HeLa cells. Through its mRNA-unwinding activity, YB-1 may thus inhibit SG assembly in cancer cells and package dormant mRNA in an unfolded state, thus preparing mRNAs for translation initiation.

Specific Recognition of a Stem-Loop RNA Structure by the Alphavirus Capsid Protein.

Alphaviruses are small enveloped viruses with positive-sense RNA genomes. During infection, the alphavirus capsid protein (Cp) selectively packages and assembles with the viral genomic RNA to form the nucleocapsid core, a process critical to the production of infectious virus. Prior studies of the alphavirus Semliki Forest virus (SFV) showed that packaging and assembly are promoted by Cp binding to multiple high affinity sites on the genomic RNA. Here, we developed an in vitro Cp binding assay based on fluorescently labeled RNA oligos. We used this assay to explore the RNA sequence and structure requirements for Cp binding to site #1, the top binding site identified on the genomic RNA during all stages of virus assembly. Our results identify a stem-loop structure that promotes specific binding of the SFV Cp to site #1 RNA. This structure is also recognized by the Cps of the related alphaviruses chikungunya virus and Ross River virus.

Hairpin Transfer-Independent Parvovirus DNA Replication Produces Infectious Virus.

Parvoviruses package a linear single-stranded DNA genome with hairpin structures at both ends. It has been thought that terminal hairpin sequences are indispensable for viral DNA replication. Here, we provide evidence that the hairpin-deleted duplex genomes of human bocavirus 1 (HBoV1) replicate in human embryonic kidney 293 (HEK293) cells. We propose an alternative model for HBoV1 DNA replication in which the leading strand can initiate strand displacement without hairpin transfer. The transfection of the HBoV1 duplex genomes that retain a minimal replication origin at the right end (OriR) but with extensive deletions in the right-end hairpin (REH) generated viruses in HEK293 cells at a level 10 to 20 times lower than that of the wild-type (WT) duplex genome. Importantly, these viruses that have a genome with various deletions after the OriR but not the one retaining only the OriR replicated in polarized human airway epithelia. We discovered that the 18-nucleotide (nt) sequence (nt 5403 to 5420) beyond the OriR was sufficient to confer virus replication in polarized human airway epithelia, although its progeny virus production was ∼5 times lower than that of the WT virus. Thus, our study demonstrates that hairpin transfer-independent productive parvovirus DNA replication can occur. IMPORTANCE Hairpin transfer-independent parvovirus replication was modeled with human bocavirus 1 (HBoV1) duplex genomes whose 5' hairpin structure was ablated by various deletions. In HEK293 cells, these duplex viral genomes with ablated 5' hairpin sequence replicated efficiently and generated viruses that productively infected polarized human airway epithelium. Thus, for the first time, we reveal a previously unknown phenomenon that productive parvovirus DNA replication does not depend on the hairpin sequence at REH to initiate rolling-hairpin DNA replication. Notably, the intermediates of viral DNA replication, as revealed by two-dimensional electrophoresis, from transfections of hairpin sequence-deleted duplex genome and full-length genome in HEK293 cells as well as from virus infection of polarized human airway epithelia are similar. Thus, the establishment of the hairpin transfer-independent parvoviral DNA replication deepens our understanding of viral DNA replication and may have implications in the development of parvovirus-based viral vectors with alternative properties.

Chloroplast genome analysis of Chrysotila dentata.

Chrysotila dentata is an ecologically important marine alga contributing to the coccolith formation. In this study, a complete chloroplast (cp DNA) genome of Chrysotila dentata was sequenced by using Illumina Hiseq and was analyzed with the help of a bioinformatics tool CPGAVAS2. The circular chloroplast genome of Chrysotila dentata has a size of 109,017 bp with two inverted repeats (IRs) regions (4513 bp each) which is a common feature in most land plants and algal species. The Chrysotila dentata cp genome consists of 61 identified protein-coding genes, 30 tRNA genes and 6 rRNAs with 21 microsatellites. The phylogenetic relationship with other select algal species revealed a close phylogeny of Chrysotila dentata with Phaeocystis antarctica. This is the first report of the cp genome analysis of genus Chrysotila and the results from this study will be helpful for understanding the genetic structure and function of chloroplast in other species of Chrysotila.

Mutation of the 5'-untranslated region stem-loop mRNA structure reduces type I collagen deposition and arterial stiffness in male obese mice.

Arterial stiffening, a characteristic feature of obesity and type 2 diabetes, contributes to the development and progression of cardiovascular diseases (CVD). Currently, no effective prophylaxis or therapeutics is available to prevent or treat arterial stiffening. A better understanding of the molecular mechanisms underlying arterial stiffening is vital to identify newer targets and strategies to reduce CVD burden. A major contributor to arterial stiffening is increased collagen deposition. In the 5'-untranslated regions of mRNAs encoding for type I collagen, an evolutionally conserved stem-loop (SL) structure plays an essential role in its stability and post-transcriptional regulation. Here, we show that feeding a high-fat/high-sucrose (HFHS) diet for 28 wk increases adiposity, insulin resistance, and blood pressure in male wild-type littermates. Moreover, arterial stiffness, assessed in vivo via aortic pulse wave velocity, and ex vivo using atomic force microscopy in aortic explants or pressure myography in isolated femoral and mesenteric arteries, was also increased in those mice. Notably, all these indices of arterial stiffness, along with collagen type I levels in the vasculature, were reduced in HFHS-fed mice harboring a mutation in the 5'SL structure, relative to wild-type littermates. This protective vascular phenotype in 5'SL-mutant mice did not associate with a reduction in insulin resistance or blood pressure. These findings implicate the 5'SL structure as a putative therapeutic target to prevent or reverse arterial stiffening and CVD associated with obesity and type 2 diabetes.NEW & NOTEWORTHY In the 5'-untranslated (UTR) regions of mRNAs encoding for type I collagen, an evolutionally conserved SL structure plays an essential role in its stability and posttranscriptional regulation. We demonstrate that a mutation of the SL mRNA structure in the 5'-UTR decreases collagen type I deposition and arterial stiffness in obese mice. Targeting this evolutionarily conserved SL structure may hold promise in the management of arterial stiffening and CVD associated with obesity and type 2 diabetes.

Data-Driven and in Silico-Assisted Design of Broad Host-Range Minimal Intrinsic Terminators Adapted for Bacteria.

Efficient transcription termination relying on intrinsic terminators is critical to maintain cell fitness by avoiding unwanted read-through in bacteria. Natural intrinsic terminator (NIT) typically appears in mRNA as a hairpin followed by approximately eight conserved uridines (U-tract) at the 3' terminus. Owing to their simple structure, small size, and protein independence, assorted NITs have been redesigned as robust tools to construct gene circuits. However, most NITs exert functions to adapt to their physiological requirements rather than the demand for building synthetic gene circuits, rendering uncertain working performance when they are constructed intact in synthetic gene circuits. Here, rather than modifying NITs, we established a data-driven and in silico-assisted (DISA) design framework to forward engineer minimal intrinsic terminators (MITs). By comprehensively analyzing 75 natural intrinsic terminators from Bacillus subtilis, we revealed that two pivotal features, the length of the U-tract and the thermodynamics of the terminator hairpin, were involved in the sequence-activity relationship (SAR) of termination efficiency (TE). As per the SAR, we leveraged DISA to fabricate an array of MITs composed of in silico-assisted designed minimal hairpins and fixed U-tracts. Most of these MITs exhibited high TE in diverse Gram-positive and Gram-negative bacteria. In contrast, the TEs of the NITs were highly varied in different hosts. Moreover, TEs of MITs were flexibly tuned over a wide range by modulating the length of the U-tract. Overall, these results demonstrate an efficient framework to forward design functional and broad host-range terminators independent of tedious and iterative screening of mutagenesis libraries of natural terminators.

Secondary Structure Libraries for Artificial Evolution Experiments.

Methods of artificial evolution such as SELEX and in vitro selection have made it possible to isolate RNA and DNA motifs with a wide range of functions from large random sequence libraries. Once the primary sequence of a functional motif is known, the sequence space around it can be comprehensively explored using a combination of random mutagenesis and selection. However, methods to explore the sequence space of a secondary structure are not as well characterized. Here we address this question by describing a method to construct libraries in a single synthesis which are enriched for sequences with the potential to form a specific secondary structure, such as that of an aptamer, ribozyme, or deoxyribozyme. Although interactions such as base pairs cannot be encoded in a library using conventional DNA synthesizers, it is possible to modulate the probability that two positions will have the potential to pair by biasing the nucleotide composition at these positions. Here we show how to maximize this probability for each of the possible ways to encode a pair (in this study defined as A-U or U-A or C-G or G-C or G.U or U.G). We then use these optimized coding schemes to calculate the number of different variants of model stems and secondary structures expected to occur in a library for a series of structures in which the number of pairs and the extent of conservation of unpaired positions is systematically varied. Our calculations reveal a tradeoff between maximizing the probability of forming a pair and maximizing the number of possible variants of a desired secondary structure that can occur in the library. They also indicate that the optimal coding strategy for a library depends on the complexity of the motif being characterized. Because this approach provides a simple way to generate libraries enriched for sequences with the potential to form a specific secondary structure, we anticipate that it should be useful for the optimization and structural characterization of functional nucleic acid motifs.

Opposing roles for Egalitarian and Staufen in transport, anchoring and localization of oskar mRNA in the Drosophila oocyte.

Localization of oskar mRNA includes two distinct phases: transport from nurse cells to the oocyte, a process typically accompanied by cortical anchoring in the oocyte, followed by posterior localization within the oocyte. Signals within the oskar 3' UTR directing transport are individually weak, a feature previously hypothesized to facilitate exchange between the different localization machineries. We show that alteration of the SL2a stem-loop structure containing the oskar transport and anchoring signal (TAS) removes an inhibitory effect such that in vitro binding by the RNA transport factor, Egalitarian, is elevated as is in vivo transport from the nurse cells into the oocyte. Cortical anchoring within the oocyte is also enhanced, interfering with posterior localization. We also show that mutation of Staufen recognized structures (SRSs), predicted binding sites for Staufen, disrupts posterior localization of oskar mRNA just as in staufen mutants. Two SRSs in SL2a, one overlapping the Egalitarian binding site, are inferred to mediate Staufen-dependent inhibition of TAS anchoring activity, thereby promoting posterior localization. The other three SRSs in the oskar 3' UTR are also required for posterior localization, including two located distant from any known transport signal. Staufen, thus, plays multiple roles in localization of oskar mRNA.

Amplification dynamics of miniature inverted-repeat transposable elements and their impact on rice trait variability.

Transposable elements (TEs) are a rich source of genetic variability. Among TEs, miniature inverted-repeat TEs (MITEs) are of particular interest as they are present in high copy numbers in plant genomes and are closely associated with genes. MITEs are deletion derivatives of class II transposons, and can be mobilized by the transposases encoded by the latter through a typical cut-and-paste mechanism. However, MITEs are typically present at much higher copy numbers than class II transposons. We present here an analysis of 103 109 transposon insertion polymorphisms (TIPs) in 738 Oryza sativa genomes representing the main rice population groups. We show that an important fraction of MITE insertions has been fixed in rice concomitantly with its domestication. However, another fraction of MITE insertions is present at low frequencies. We performed MITE TIP-genome-wide association studies (TIP-GWAS) to study the impact of these elements on agronomically important traits and found that these elements uncover more trait associations than single nucleotide polymorphisms (SNPs) on important phenotypes such as grain width. Finally, using SNP-GWAS and TIP-GWAS we provide evidence of the replicative amplification of MITEs.

1H, 13C, 15N and 31P chemical shift assignment for stem-loop 4 from the 5'-UTR of SARS-CoV-2.

The SARS-CoV-2 virus is the cause of the respiratory disease COVID-19. As of today, therapeutic interventions in severe COVID-19 cases are still not available as no effective therapeutics have been developed so far. Despite the ongoing development of a number of effective vaccines, therapeutics to fight the disease once it has been contracted will still be required. Promising targets for the development of antiviral agents against SARS-CoV-2 can be found in the viral RNA genome. The 5'- and 3'-genomic ends of the 30 kb SCoV-2 genome are highly conserved among Betacoronaviruses and contain structured RNA elements involved in the translation and replication of the viral genome. The 40 nucleotides (nt) long highly conserved stem-loop 4 (5_SL4) is located within the 5'-untranslated region (5'-UTR) important for viral replication. 5_SL4 features an extended stem structure disrupted by several pyrimidine mismatches and is capped by a pentaloop. Here, we report extensive 1H, 13C, 15N and 31P resonance assignments of 5_SL4 as the basis for in-depth structural and ligand screening studies by solution NMR spectroscopy.

Live imaging and quantitation of nascent transcription using the MS2/MCP system in the Drosophila embryo.

Visualizing transcription live in Drosophila is providing important new insights into the spatiotemporal regulation of transcription. Here, we describe a protocol to visualize and quantitate transcription from gene loci that are tagged with MS2 stem-loop sequences in the Drosophila embryo. MS2 stem-loop sequences are recognized by a coat protein fused to a fluorescent protein and visualized with microscopy. We also describe an analysis pipeline to extract and subsequently quantify transcription dynamics. For complete details on the use and execution of this protocol, please refer to Hoppe et al. (2020).

CRISPR-Cas9 strategies to insert MS2 stem-loops into endogenous loci in Drosophila embryos.

CRISPR-Cas9 genome editing has transformed biology by enabling site-specific genome modifications to be simply engineered. Here, we describe two CRISPR-Cas9 approaches to introduce MS2 stem-loop sequences into endogenous gene loci in Drosophila. This can facilitate live imaging of nascent transcription in Drosophila. For complete details on the use and execution of this protocol, please refer to Hoppe et al. (2020).

Genome-wide analysis of long noncoding RNAs, 24-nt siRNAs, DNA methylation and H3K27me3 marks in Brassica rapa.

Long noncoding RNAs (lncRNAs) are RNA fragments that generally do not code for a protein but are involved in epigenetic gene regulation. In this study, lncRNAs of Brassica rapa were classified into long intergenic noncoding RNAs, natural antisense RNAs, and intronic noncoding RNAs and their expression analyzed in relation to genome-wide 24-nt small interfering RNAs (siRNAs), DNA methylation, and histone H3 lysine 27 trimethylation marks (H3K27me3). More than 65% of the lncRNAs analyzed consisted of one exon, and more than 55% overlapped with inverted repeat regions (IRRs). Overlap of lncRNAs with IRRs or genomic regions encoding for 24-nt siRNAs resulted in increased DNA methylation levels when both were present. LncRNA did not overlap greatly with H3K27me3 marks, but the expression level of intronic noncoding RNAs that did coincide with H3K27me3 marks was higher than without H3K27me3 marks. The Brassica genus comprises important vegetables and oil seed crops grown across the world. B. rapa is a diploid (AA genome) thought to be one of the ancestral species of both B. juncea (AABB genome) and B. napus (AACC) through genome merging (allotetrapolyploidization). Complex genome restructuring and epigenetic alterations are thought to be involved in these allotetrapolyploidization events. Comparison of lncRNAs between B. rapa and B. nigra, B. oleracea, B. juncea, and B. napus showed the highest conservation with B. oleracea. This study presents a comprehensive analysis of the epigenome structure of B. rapa at multi-epigenetic levels (siRNAs, DNA methylation, H3K27me3, and lncRNAs) and identified a suite of candidate lncRNAs that may be epigenetically regulated in the Brassica genus.

Uncovering dynamic evolution in the plastid genome of seven Ligusticum species provides insights into species discrimination and phylogenetic implications.

Ligusticum L., one of the largest members in Apiaceae, encompasses medicinally important plants, the taxonomic statuses of which have been proved to be difficult to resolve. In the current study, the complete chloroplast genomes of seven crucial plants of the best-known herbs in Ligusticum were presented. The seven genomes ranged from 148,275 to 148,564 bp in length with a highly conserved gene content, gene order and genomic arrangement. A shared dramatic decrease in genome size resulted from a lineage-specific inverted repeat (IR) contraction, which could potentially be a promising diagnostic character for taxonomic investigation of Ligusticum, was discovered, without affecting the synonymous rate. Although a higher variability was uncovered in hotspot divergence regions that were unevenly distributed across the chloroplast genome, a concatenated strategy for rapid species identification was proposed because separate fragments inadequately provided variation for fine resolution. Phylogenetic inference using plastid genome-scale data produced a concordant topology receiving a robust support value, which revealed that L. chuanxiong had a closer relationship with L. jeholense than L. sinense, and L. sinense cv. Fuxiong had a closer relationship to L. sinense than L. chuanxiong, for the first time. Our results not only furnish concrete evidence for clarifying Ligusticum taxonomy but also provide a solid foundation for further pharmaphylogenetic investigation.

Hairpin-Bisulfite PCR.

Ligation of a hairpin oligonucleotide to genomic DNA prior to bisulfite conversion and PCR amplification physically links the two complementary DNA strands. This additional step in the conversion procedure overcomes the limitations of conventional bisulfite sequencing where information of the cytosine methylation status is only obtained from one of the two strands of an individual DNA molecule. Sequences derived from hairpin bisulfite PCR products reveal the dynamics of this epigenetic memory system on both strands of individual DNA molecules. The chapter describes a reliable step-by-step procedure to generate hairpin-linked DNA. It also provides a guide for efficient bisulfite conversion that is suitable for both conventional and hairpin bisulfite sequencing approaches.